Single- and Multiple-Shell Uniform Sampling Schemes for Diffusion MRI Using Spherical Codes

In diffusion MRI (dMRI), a good sampling scheme is important for efficient acquisition and robust reconstruction. Diffusion weighted signal is normally acquired on single or multiple shells in q-space. Signal samples are typically distributed uniformly on different shells to make them invariant to the orientation of structures within tissue, or the laboratory coordinate frame. The Electrostatic Energy Minimization (EEM) method, originally proposed for single shell sampling scheme in dMRI, was recently generalized to multi-shell schemes, called Generalized EEM (GEEM). GEEM has been successfully used in the Human Connectome Project (HCP). However, EEM does not directly address the goal of optimal sampling, i.e., achieving large angular separation between sampling points. In this paper, we propose a more natural formulation, called Spherical Code (SC), to directly maximize the minimal angle between different samples in single or multiple shells. We consider not only continuous problems to design single or multiple shell sampling schemes, but also discrete problems to uniformly extract sub-sampled schemes from an existing single or multiple shell scheme, and to order samples in an existing scheme. We propose five algorithms to solve the above problems, including an incremental SC (ISC), a sophisticated greedy algorithm called Iterative Maximum Overlap Construction (IMOC), an 1-Opt greedy method, a Mixed Integer Linear Programming (MILP) method, and a Constrained Non-Linear Optimization (CNLO) method. To our knowledge, this is the first work to use the SC formulation for single or multiple shell sampling schemes in dMRI. Experimental results indicate that SC methods obtain larger angular separation and better rotational invariance than the state-of-the-art EEM and GEEM. The related codes and a tutorial have been released in DMRITool.Comment: Accepted by IEEE transactions on Medical Imaging. Codes have been released in dmritool https://diffusionmritool.github.io/tutorial_qspacesampling.htm

Exploring Local White Matter Geometric Structure in diffusion MRI Using Director Field Analysis

International audienceIn this abstract, inspired by microscopic theoretical treatment of phases in liquid crystals, we introduce a novel mathematical framework, called Director Field Analysis (DFA), to study local geometric structural information of white matter. DFA extracts some meaningful scalar indices related with the degree of orientational alignment, dispersion, and orientational distortion, from the Orientation Distribution Function (ODF) field reconstructed by Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI)

Director Field Analysis to Explore Local White Matter Geometric Structure in diffusion MRI

International audienceIn diffusion MRI, a tensor field or a spherical function field, e.g., an Orientation Distribution Function (ODF) field, are estimated from measured diffusion weighted images. In this paper, inspired by microscopic theoretical treatment of phases in liquid crystals, we introduce a novel mathematical framework, called Director Field Analysis (DFA), to study local geometric structural information of white matter from the estimated tensor field or spherical function field. 1) We propose Orientational Order (OO) and Orientational Dispersion (OD) indices to describe the degree of alignment and dispersion of a spherical function in each voxel; 2) We estimate a local orthogonal coordinate frame in each voxel with anisotropic diffusion; 3) Finally, we define three indices to describe three types of orientational distortion (splay, bend, and twist) in a local spatial neighborhood, and a total distortion index to describe distortions of all three types. To our knowledge, this is the first work to quantitatively describe orientational distortion (splay, bend, and twist) in diffusion MRI. The proposed scalar indices are useful to detect local geometric changes of white matter for voxel-based or tract-based analysis in both DTI and HARDI acquisitions

High Resolution Ex Vivo Diffusion Tensor Distribution MRI of Neural Tissue

Neural tissue microstructure plays a key role in developmental, physiological and pathophysiological processes. In the continuing quest to characterize it at ever finer length scales, we use a novel diffusion tensor distribution (DTD) paradigm to probe microstructural features much smaller than the nominal MRI voxel size. We first assume the DTD is a normal tensor variate distribution constrained to lie on the manifold of positive definite matrices, characterized by a mean and covariance tensor. We then estimate the DTD using Monte Carlo signal inversion combined with parsimonious model selection framework that exploits a hierarchy of symmetries of mean and covariance tensors. High resolution multiple pulsed field gradient (mPFG) MRI measurements were performed on a homogeneous isotropic diffusion phantom (PDMS) for control, and excised visual cortex and spinal cord of macaque monkey to investigate the capabilities of DTD MRI in revealing neural tissue microstructural features using strong gradients not typically available in clinical MRI scanners. DTD-derived stains and glyphs, which disentangle size, shape, and orientation heterogeneities of microscopic diffusion tensors, are presented for all samples along with the distribution of the mean diffusivity (MD) within each voxel. We also present a new glyph to visualize the symmetric (kurtosis) and asymmetric parts of the fourth-order covariance tensor. An isotropic mean diffusion tensor and zero covariance tensor was found for the isotropic PDMS phantom, as expected, while the covariance tensor (both symmetric and asymmetric parts) for neural tissue was non-zero indicating that the kurtosis tensor may not be sufficient to fully describe the microstructure. Cortical layers were clearly delineated in the higher moments of the MD spectrum consistent with histology, and microscopic anisotropy was detected in both gray and white matter of neural tissue. DTD MRI captures crossing and splaying white matter fibers penetrating into the cortex, and skewed fiber diameter distributions in the white matter tracts within the cortex and spinal cord. DTD MRI was also shown to subsume diffusion tensor imaging (DTI) while providing additional microstructural information about tissue heterogeneity and microscopic anisotropy within each voxel.Peer reviewe

Multi-site genetic analysis of diffusion images and voxelwise heritability analysis : a pilot project of the ENIGMA–DTI working group

The ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) Consortium was set up to analyze brain measures and genotypes from multiple sites across the world to improve the power to detect genetic variants that influence the brain. Diffusion tensor imaging (DTI) yields quantitative measures sensitive to brain development and degeneration, and some common genetic variants may be associated with white matter integrity or connectivity. DTI measures, such as the fractional anisotropy (FA) of water diffusion, may be useful for identifying genetic variants that influence brain microstructure. However, genome-wide association studies (GWAS) require large populations to obtain sufficient power to detect and replicate significant effects, motivating a multi-site consortium effort. As part of an ENIGMA–DTI working group, we analyzed high-resolution FA images from multiple imaging sites across North America, Australia, and Europe, to address the challenge of harmonizing imaging data collected at multiple sites. Four hundred images of healthy adults aged 18–85 from four sites were used to create a template and corresponding skeletonized FA image as a common reference space. Using twin and pedigree samples of different ethnicities, we used our common template to evaluate the heritability of tract-derived FA measures. We show that our template is reliable for integrating multiple datasets by combining results through meta-analysis and unifying the data through exploratory mega-analyses. Our results may help prioritize regions of the FA map that are consistently influenced by additive genetic factors for future genetic discovery studies. Protocols and templates are publicly available at (http://enigma.loni.ucla.edu/ongoing/dti-working-group/)

Heterogeneity of fractional anisotropy and mean diffusivity measurements by in vivo diffusion tensor imaging in normal human hearts

Background: Cardiac diffusion tensor imaging (cDTI) by cardiovascular magnetic resonance has the potential to assess microstructural changes through measures of fractional anisotropy (FA) and mean diffusivity (MD). However, normal variation in regional and transmural FA and MD is not well described. Methods: Twenty normal subjects were scanned using an optimised cDTI sequence at 3T in systole. FA and MD were quantified in 3 transmural layers and 4 regional myocardial walls. Results: FA was higher in the mesocardium (0.46 ±0.04) than the endocardium (0.40 ±0.04, p≤0.001) and epicardium (0.39 ±0.04, p≤0.001). On regional analysis, the FA in the septum was greater than the lateral wall (0.44 ±0.03 vs 0.40 ±0.05 p = 0.04). There was a transmural gradient in MD increasing towards the endocardium (epicardium 0.87 ±0.07 vs endocardium 0.91 ±0.08×10-3 mm2/s, p = 0.04). With the lateral wall (0.87 ± 0.08×10-3 mm2/s) as the reference, the MD was higher in the anterior wall (0.92 ±0.08×10-3 mm2/s, p = 0.016) and septum (0.92 ±0.07×10-3 mm2/s, p = 0.028). Transmurally the signal to noise ratio (SNR) was greatest in the mesocardium (14.5 ±2.5 vs endocardium 13.1 ±2.2, p<0.001; vs epicardium 12.0 ± 2.4, p<0.001) and regionally in the septum (16.0 ±3.4 vs lateral wall 11.5 ± 1.5, p<0.001). Transmural analysis suggested a relative reduction in the rate of change in helical angle (HA) within the mesocardium. Conclusions: In vivo FA and MD measurements in normal human heart are heterogeneous, varying significantly transmurally and regionally. Contributors to this heterogeneity are many, complex and interactive, but include SNR, variations in cardiac microstructure, partial volume effects and strain. These data indicate that the potential clinical use of FA and MD would require measurement standardisation by myocardial region and layer, unless pathological changes substantially exceed the normal variation identified